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In isotropic systems below the transition temperature, the massless Goldstone 
modes imply critical infrared singularities in the statics and dynamics along the en- 
tire coexistence curve. We examine the important question whether these coexistence 
anomalies are of relevance also in more realistic systems displaying anisotropies. By ap- 
plying a generalized renormalization scheme to the time-dependent Ginzburg-Landau 
models, we treat two quite different but characteristic cases, namely the influence of 
(i) weak cubic anisotropies, and (ii) long-range dipolar interactions. In the presence of 
cubic terms, the transverse excitations acquire a mass and thus one expects the theory 
to approach an uncritical "Gaussian" regime in the limit q — > and uj — > 0. Therefore, 
we first consider the one-component case in order to show that our formalism also pro- 
vides a consistent description of the crossover into an asymptotically Gaussian theory. 
In the case of (weak) cubic anisotropies, the fact that for n < 4 the fluctuations tend 
to restore the 0(n)-symmetry at the critical point proves to be most important, since 
under these circumstances coexistence-type singularities may be found in an interme- 
diate wavenumber and frequency range. The dipolar interaction induces an anisotropy 
in momentum space, which does not completely destroy the massless character of the 
transverse fluctuations, but only effectively reduces the number of Goldstone modes by 
one. Remarkably, similar to the isotropic case the asymptotic theory can be treated 
exactly. For n > 3 we find coexistence anomalies governed by the isotropic power laws. 
However, the amplitudes of the respective scaling functions depend on the angle between 
order parameter and external wavevector. 
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I. INTRODUCTION 



In this paper, we study the static and dynamic susceptibilities of the n-component 
time-dependent Ginzburg-Landau models (also called relaxational models A and B ac- 
cording to Hohenberg and Halperin's classification [1]) below the transition temperature, 
where in addition to an underlying 0(n)-symmetry of the order parameter field two char- 
acteristic anisotropies are taken into account, namely (i) weak cubic anharmonic terms, 
and (ii) long-range dipolar interactions. 

Our considerations have the following background. In the ordered phase of ide- 
ally isotropic systems, as a consequence of the spontaneously broken symmetry there 
appear n — 1 massless Goldstone modes which lead to infrared singularities in certain 
correlation functions for all temperatures below the critical temperature T c . These so- 
called coexistence anomalies have been studied on the basis of several exponentiation 
and renormalization schemes for the static </> 4 -model [2-9] as well as the dynamic re- 
laxational models [5,10]. Lawrie's approach [8] to take into account both critical and 
Goldstone fluctuations within a suitable crossover theory, has recently been extended by 
the present authors to the 0(n) -symmetric time-dependent Ginzburg-Landau models 
[11]. Generally, our calculations are based on the path-integral formulation of dynami- 
cal perturbation theory for Langevin-type equations of motion as developed by Bausch, 
Janssen, Wagner [12] and De Dominicis [13], supplemented by Amit and Goldschmidt's 
generalized renormalization procedure especially designed for crossover phenomena [14]. 
To our opinion, the formalism is rather transparent, the exact results in the coexistence 
limit can be easily implemented, and the connections to different approaches, e.g. the 
1 /n-expansion [10,15,16], can be very clearly demonstrated. One of the major addi- 
tional advantages is the possible applicability to even more complicated models. Some 
of the central results for the dynamic susceptibilities and correlation functions of the 
relaxational models A and B are listed in Table I. We also remark that lately a direct 
comparison with experimental findings of ultrasonic attenuation measurements within 
the incommensurate phase of the A 2 -BX 4 -family crystals has been accomplished [17]. 

Having established a theory for ideally 0(n)-symmetric systems, it is quite natural 
a question to ask whether the striking infrared anomalies found in the coexistence regime 
completely disappear when (weak) anisotropies are introduced, or certain remnants of 
the Goldstone excitations may still be experimentally relevant. This issue is even more 
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important as for many cases, especially for structural and magnetic phase transforma- 
tions with an underlying discrete lattice symmetry, the assumption of an isotropic order 
parameter field is certainly an idealization which is only approximately justified. On 
the other hand, one would expect a smooth transition of the physical properties from 
an (9(n)-symmetric model system to a real, anisotropic crystal. Wallace [18] and Nelson 
[4] have noted that below but close to T c even in weakly anisotropic sytems very strong 
fluctuations induced by quasi- Goldstone modes should be prominent. The dynamics in 
the high-frequency region of such cubic crystals in the ordered phase was studied by 
Meissner, Menyhard, and Szepfalusy [19]. Our focus, however, will be the crossover 
from static and dynamic critical behavior into the hydrodynamic regime of vanishing 
external wavenumbers and frequencies. 

The outline of our paper is as follows. In the following Section we briefly review 
the main results obtained for the (9(n)-symmetric relaxational models below the tran- 
sition temperature. The information presented here only serves for comparison with 
the chapters below. For a thorough explanation and discussion of the renormalization 
technique used, we refer the interested reader to Ref.fll]. By specializing to the case 
n = 1, we shall demonstrate in Section III that our procedure also provides a consistent 
description of a crossover to an asymptotically uncritical theory which can be described 
by a Gaussian model. This constitutes an inevitable preliminary, as the existence of 
cubic anisotropies in the free energy for the order parameter field will eventually (i.e.: 
for q — > and to — > 0) cause all the fluctuations to cease. Under certain circum- 
stances, however, namely if in the vicinity of the critical point the 0(n) -symmetry is 
dynamically restored, coexistence-type singularities may be found in an intermediate 
wavenumber and frequency region, as we are going to point out in Section IV. There, 
for n < 4 the crossover from the critical behavior to the asymptotically Gaussian model 
via a Goldstone-mode dominated regime will be discussed in detail, based on a one- 
loop approximation. In Section V we shall turn to our second generic example, namely 
long-range dipolar interactions introducing anisotropies in momentum space. Our inves- 
tigation of the asymptotic theory, which remarkably can be treated exactly on very much 
the same basis as the isotropic model, shows that effectively the number of massless ex- 
citations is reduced by one. Hence for n > 3 coexistence anomalies persist, characterized 
by just the same power laws as summarized in Table I. However, the amplitudes of the 
associated scaling functions are generally lower, the reduction factors depending on the 
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angle between the direction of the static order parameter and the external wavevector. 
At last we shall discuss and summarize our results. In the Appendix we list the dia- 
grams and the corresponding analytical results for the two-point vertex functions up 
to one-loop order for the time-dependent Ginzburg-Landau models taking into account 
cubic terms. 



II. COEXISTENCE ANOMALIES IN THE ISOTROPIC 
RELAXATIONAL MODELS 

For convenience we start with a brief review of the static and dynamic properties 
of the 0(n)-symmetric time-dependent Ginzburg-Landau models below the transition 
temperature. We shall merely repeat the basic definitions of the models under con- 
sideration and present the central results of the theory for later comparison with the 
anisotropic cases. Further details, especially concerning the dynamical perturbation 
expansion in the ordered phase and the generalized renormalization scheme used to 
describe the crossover behavior, may be found in our recent paper [11]. 

A. Model equations and perturbation theory 

We are interested in the statics and dynamics near a second-order phase transition 
which phenomenologically can be described by an n-component order parameter field 
4>q , a = 1, . . . , n, and the following 0(n) -symmetric expansion of the Ginzburg-Landau 
free-energy functional ("Hamiltonian") with respect to </>q, 

2" 

ct=l a=l \a=l 

The parameter tq is proportional to the separation from the mean-field transition tem- 
perature T®, and hence can be written as a sum of the fluctuation-induced T c -shift r$ c 
and the reduced temperature variable with respect to the critical point T c 

r = r 0c + ^r^ . (2.16) 

The positive coupling uq defines the strength of the isotropic anharmonic term. In 
thermal equilibrium, the probability density for a specific configuration <j>o(x) is 

P ^o}] = j v[m]e - Hm}] ■ ( 2 - 2 ) 
4 
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To constitute the dynamics of our model we assume that (i) the only "slow" variable 
in the system is the order parameter itself, i.e.: there is no conserved quantity besides 
perhaps (f>Q, and (ii) no reversible forces appear in the generalized Langevin equation, 
which is the simplest conceivable case. Therefore we have purely relaxational behavior 
where the damping oc Ao may be considered as originating from the action of the " fast" 
degrees of freedom subsummed in the fluctuating forces ( a . Thus our basic equation of 
motion reads 

|«( 2 , ( ) = - M ^.*ggH + <-(*,«) . (2.3) 

Here two situations have to be distinguished, namely either there is no conservation law 
for the order parameter field, and the system simply relaxes into its equlibrium state 
after a distortion [a = in Eq.(2.3)], or 4>q is a conserved quantity itself and hence 
follows a diffusion equation (a = 2). These two cases are referred to as model A or 
model B, respectively, according to the classification by Hohenberg and Halperin [1]. 

If one takes the probability distribution of the stochastic forces ( a to be of a 
Gaussian form, one finds for the first two moments 

(C(x,t)) = (2.4a) 
{C(x,t)C P (x,t')} = 2X (tV) a S ap 5(x-x')S(t-t') , (2.46) 

which completes the definition of the time-dependent Ginzburg-Landau models. The 
Einstein relation (2.4b) ensures that the equilibrium probability density for a certain 
configuration 4>q(x) is indeed given by Eq.(2.2). 

Dynamical correlation functions are now calculated within the field-theoretical ap- 
proach to generalized Langevin equations as developed by Bausch, Janssen, Wagner 
[12] and De Dominicis [13]. Starting from the assumed Gaussian distribution for the 
stochastic forces, one eliminates the ( a using the equation of motion (2.3). After a Gaus- 
sian transformation and the accompanying introduction of Martin-Siggia-Rose auxiliary 
fields 0o , the probability density for the order parameter fluctuations 0q finally reads 

p[{<t>d }] = J mk}]p[{k},{m « J v[{i^}\e j ^^ , (2 . 5 a) 

where the statistical weight for a specific configuration is determined by the so-called 
Janssen-De Dominicis functional J. For the time-dependent Ginzburg-Landau models 
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one finds 

Jim, m] = Jd d xJdtj2 [* A (• vr * - * + A (i vr 5H[ ^I }] ) 

(2.56) 

here we have omitted a term stemming from the functional derivative of the transfor- 
mation from the ( a to the 0q, because it exactly cancels those contributions to the 
perturbation series which violate causality [12]. 

Based on the path-integral formulation (2.5) for non-linear Langevin dynamics, a 
perturbation expansion with respect to the anharmonicity oc uq and its diagrammatic 
representation by propagators and vertices can be constructed in complete analogy to 
the static case (see e.g. Ref.[20]). As usual, one defines a generating functional for Ap- 
point Green functions and cumulants by inserting source terms into (2.5b); furthermore, 
after a Legendre transform one may also obtain the vertex functions by appropriate 
functional derivatives. Details may be found in Refs. [20,12,11]. 

In order to evaluate the dynamic (linear) response functions within our formalism, 
we have to take into account the action of an external field h a . This leads to an 
additional term — J2 a h a 4>q in the Hamiltonian (2.1a), and thus to a modified Janssen- 
De Dominicis functional [12] 

A{fch{<P5}] = J[{4tih {<%}]+ fd d x fdtJ2~h a ^(i^) a ^ ■ (2-6) 

J J a 

Hence the dynamic susceptibility is given by the expression 



a/3,* + . -/ ^ _ *(#)0M)) 



an / -> , -,/ ,/\ 



= \ Q (cj>%{x,t){iV) a ~4{x',t')) , (2.7a) 



6hP(x',t') 

and its Fourier transform is intimately connected with the so-called response propagator 

X f(q,u) = X q a G ^ (j)f} (q,oj) . (2.7b) 

The dynamical correlation function can either be calculated directly or derived from the 
following (classical) fluctuation-dissipation theorem [12] 

2 Ao q a 

G 'o4> a 4>p (q, u) = — ImG i a ^(g» . (2.8) 
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B. Renormalization in the ordered phase and coexistence anomalies 

Below the critical temperature T c a spontaneous order parameter 0o appears, which 
we assume to point in the n-th direction of the order parameter space. For convenience 
we define new transverse (a = 1, . . . ,n — 1) and longitudinal fields 



a oJ \<PoJ V°"o + </>o 



(2.9a) 



with vanishing thermal average 

«) = (a ) = . (2.96) 

In addition we shall use the parametrization 
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mo , (2.10) 

where the quantity mo has the dimension of a mass or inverse length. By evaluating 
the condition (2.9b) for the longitudinal field one arrives at the equation of state 

r + ^=A + O(u 2 ) , (2.11) 

which allows the elimination of the temperature variable ro in favor of the parameter 
mo [11,8]; to express the perturbation expansion in terms of mo has the considerable 
advantage that the T c -shift is already properly included (compare the discussion in 
Ref.[21]). 

From a general dynamic Ward-Takahashi identity [11] the following exact relation 
for the transverse two-point vertex function can be derived 

$ o T o * ir (q = 0,u = 0) = h n , (2.12) 

implying the massless character of the n — 1 transverse excitations for vanishing external 
field h n — > 0. These Goldstone modes induce infrared singularities that are not restricted 
to the vicinity of the critical point, but persist in the entire ordered phase. 

Very remarkably, however, in contrast to the critical region (mo ~ 0) the coexistence 
regime T < T c , q — > 0, u> — > 0, h n — > allows for an exact treatment [8,15,11]. 
Heuristically, the above limits may be replaced by fixed wavenumbers and frequencies 
and mo — > oo instead. The renormalization group analysis a posteriori confirms the 



assertion that the infrared behavior is indeed characterized by a diverging longitudinal 
mass parameter m. By introducing new longitudinal fields 



Q = l 



<Po{q,u) = m a (q,uj)+ / V ttq (<f, u/) ttq (g - <f, u -u) + J— A 5(q) 5(u) 

b Ja'Ju' ~t V «0 

(2.13a) 



o (£ O j) = m o c- o (q 1 u) + -^ f [ ^Y^^tf^') Kit (2-136) 



5 a=l 

(provided q ^ in the case of model B) and performing the limit mo — > oo — under 
the assumption uq/tuq ^0 — one arrives at a purely harmonic asymptotic functional. 
(Here we have introduced the short-hand notation J q J u ■ • • = ^ 2 J) d +- 
In the coexistence regime merely transverse fluctuations exist, and the corresponding 
Goldstone propagators are identical to those from mean-field theory, as has been an- 
ticipated in earlier work [2,3,5]. On the other hand, by using (2.13) one finds for the 
(original) longitudinal correlation functions in the asymptotic limit [11] exact results 
(see Eq.(3.18) of Ref.fll]) that correspond to (i) the one- loop expressions for the two- 
point cumulants, and (ii) the leading order of the 1 /n-expansion for the two-point vertex 
functions [15,16], which can be summed up via a geometric series (see Fig. 2 of Ref.fll]). 
In the asymptotic model the ultraviolet divergences can be removed by introducing 
renormalized counterparts m and u for the order parameter and coupling constant, re- 
spectively, where m oc u with a finite proportionality factor. Within the renormalized 
theory, the remaining infrared divergences are transcribed to an anomalous dimension 
of the mass m which eventually yields the correct behavior for low wavenumbers and 
frequencies. 

Following Lawrie's work [8], we apply Amit and Goldschmidt's generalized minimal 
subtraction procedure [14] in order to describe the entire crossover from the critical 
behavior to the Goldstone regime. However, as we do not want to weaken the exact 
statements valid in the coexistence limit by introducing any further approximations, we 
refrain from e-expansion [11], where 

0< e = 4-d<2 (2.14) 

is the difference to the upper critical dimension of the -model. Following the argu- 
ments of Schloms and Dohm [21], quite generally the extraction of the correct infrared 
behavior via the analysis of the ultraviolet divergences appearing as e-poles within the 



dimensional regularization scheme does not require an expansion with respect to the 
distance from the upper critical dimension. The latter is needed if one wants to de- 
fine a "small" expansion parameter in the perturbation series which can be avoided by 
a suitable resummation procedure. Our case is even simpler, though, for the theory 
asymptotically reduces to the one-loop expressions. 

Neglecting the Fisher corrections connected with the static exponent rj we thus 
use a one-loop approximation for the entire wavenumber and frequency region. On 
this level we can drop any field renormalizations; furthermore, because the Z-factor for 
the Onsager coefficient as a consequence of the fluctuation-dissipation theorem (2.8) 
can be expressed in terms of and [11,12] (see also Eq.(5.18) below), one finds 
Z\ = 1 + 0(ul ) for the relaxational models. Hence the only renormalized quantities we 
have to introduce are 



m 



u 



r?-\ 2 -2 

Z m m v 
Z~ x u Q A d \i~ 



(2.15a) 
(2.156) 



where for convenience we have explicitly separated the naive dimensions and the geo- 
metric factor [21] 



A d = S d T 3 



d 



d 



1 = 



r(3-|) 



(2.16) 



2 ) 2 d - 2 ni (d-2) 
From the explicit one-loop results for the Z-factors defined in (2.15) we proceed to 
calculate Wilson's flow functions 



( m (u, m) = /I 
p u (u,m) = [i 



d_ 
d_ 



In 



u 



m 



rat 



-2 + 



n 
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1 3 
-u + 



u 



u 



n-1 
— e H — u + 



2 (l + m 2 )^/ 2 
3 u 



(2.17a) 
(2.176) 



6 " ' 2 (l + m 2 ) 1 +"/\ 
They enter the renormalization group equations for the cumulants, which are solved via 
the introduction of the characteristics fj>(£) = \it with the result 



G% N <ji, A, m, u, M, M) = G%M A, m(£), u(£), {q/^i}, {u/{^) 2 }) 



(2.18) 



where the flow- dependent parameters are to be determined as solutions of the coupled 
first-order differential equations 



t 



dm{t) 
d£ 



(2.19a) 
(2.196) 



with the initial conditions m(l) = m, u(l) = u. 

The fixed points of the renormalization group characterizing the asymptotic behav- 
ior are given by the zeros of the /9-function (2.19b). At the critical point (m = 0) the 
Heisenberg fixed point 

u* H = — — (2.20a) 
H n+8 V ; 

becomes stable, and the associated anomalous dimension for the mass parameter is 

C mH = -2 + e . (2.206) 

On the contrary, in the limit m — > oo Lawrie's coexistence fixed point [8] is approached 

u% = — — (2.21a) 
c n-1 V ' 

C mC = -2 + e ■ (2.216) 

hence indeed m(£) 2 oc £~ 2+e — > oo for £ — > in the range < e < 2, and the above 
treatment of the asymptotic model becomes meaningful. For d > 4, the Gaussian fixed 
point is stable and the coexistence anomalies (as well as any non-trivial critical behavior) 
disappear. At the critical dimension d c = 4 there will be logarithmic corrections. On 
approaching d — > 2, m(£) no longer diverges and the assumption of a homogeneous long- 
range order parameter turns out to be inconsistent, which is to be expected according 
to the theorem by Mermin, Wagner [22], and Hohenberg [23]. 

In Ref.fll] approximate analytical as well as numerical solutions of the flow equa- 
tions (2.19) are discussed in detail (see Fig. 5 for y = 0). Here we just quote that because 
of Z m — Z u the ratio 

m{£? g2 _ £ _ m(l) 2 

u{£) 1 ~ nil) [ ] 

is invariant under the renormalization group, and the important property that for small 
initial values m(l) < 0.1, the flows coincide when expressed in terms of the scaling 
variable 

X = m(l) 2 /( 2 -) ' (2 ' 23) 

which turns out to be the origin for the general scaling laws obeyed by the static and 
dynamic susceptibilities. In order to determine the scaling functions we finally have to 
use a matching condition fixing at least one of the arguments on the right-hand side of 
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(2.18b) at a sufficiently large value such that a conventional perturbation approach is 
applicable. In this context a very suitable choice will be 

,2 



e = 



(2.24) 



i.e.: £ 2 | x^ 1 (q, u) | is precisely identified with the distance to the critical surface of 
the Goldstone singularities. 

Taking advantage of the scaling property (2.23) and inserting (2.24) into the per- 
turbational expressions for the dynamic susceptibilities, one finds that the response 
functions are subject to the general scaling law 

x -) L {r,q,u) = f-VxvM'U/uc) , (2-25) 

where 

r 1 = A*™(l) 2/(2 - e+J?) oc| t r (2.26a) 

is the inverse static correlation length [remember that m(l) oc oc| r \@ and the scaling 
relation v = 2 (3/(2 — e + rj)] with rj = to one-loop order, and 

uj c = \ix 2 ~ z i z oc| r \ zu (2.266) 

denotes the characteristic frequency scale. The asymptotic power laws for the coexis- 
tence singularities are fixed by the anomalous dimension (2.21b); remarkably one finds 
to leading order for q — > and u> — > 

XL (q,oj) oc 2 \xAq^)\ e/2 • (2.27) 

This means that the longitudinal correlation function displays an anomalous line-shape, 
considerably stretched on the frequency scale in comparison with a Lorentzian, and 
corresponding to an algebraic decay of GL(x,t) for long times. Of course, Eq.(2.24) can 
be specialized to the static limit (ui = 0) or the case q^0,uj/q a = const, for which the 
resulting coexistence anomalies are summarized in Table I. For more information about 
the crossover features compare Figs. 7-9 and the discussion in Ref.fll]. 
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III. CROSSOVER TO AN ASYMPTOTICALLY GAUSSIAN 
THEORY: THE CASE n = 1 



In the case of a single-component order parameter there will be no massless modes 
in the low-temperature phase, of course, and hence infrared singularities are confined 
to the vicinity of T c . However, one expects not only the critical region itself, but also 
the crossover into the asymptotic "Gaussian" regime to display a universal character. 
As an example, we therefore extrapolate our generalized renormalization scheme to 
the situation T < T c for n = 1, because the formalism allows the separation of the 
universal behavior from any features depending on the specific initial values of the flow 
parameters. This approach to the crossover towards a Gaussian model will turn out to 
be a necessary prerequisite for an appropriate discussion of the anisotropic n-component 
relaxational models. 

Hence we put n = 1 in the flow equations (2.19) using the explicit one-loop results 
(2.17). Again, the m-dependent term describes the freezing-out of the fluctuations on 
leaving the critical region. For m — > oo only the mean-field contributions to Wilson's 
(- and /^-functions remain, leading to a power-law behavior of the flow parameters 
according to their canonical dimensions 

m^aT 1 (3.1a) 
u(£) oc £~ e . (3.16) 

Thus the anharmonicity u(£) diverges (see also Eq. (2.21a) for n — > 1); however, for large 
m{£) the relevant effective coupling entering the flow functions vanishes according to 

This eventually leads to mean-field exponents, which is precisely what should be ex- 
pected for a Gaussian model. On the other hand, the loop-corrections to the scaling 
functions are of the order um~ e oc £°, i.e.: a non-universal constant. Hence the asymp- 
totic perturbation theory, although non-trivial with respect to the anharmonicity, does 
not contain any infrared singularities. Altogether, this seems to provide a consistent 
description for the crossover to an asymptotically uncritical regime. 

These considerations can be confirmed by the investigation of an approximate so- 
lution of Eq.(2.19b) simply using (3.1a) for the complete crossover region. Instead of 
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Eqs.(4.26a), (4.27a) of Ref.[ll] one finds for e > 1 

u£~ e 

U ^ = 1 + g [( m 2 + £2)- e /2 _ (m 2 + ' ^ 3a ) 

while in the case e = the result is identical to Eq. (4.27b) of Ref.fll] for n = 1 



u 

' m 2 + l 



Hence at d c = 4 even the logarithmic corrections in the ordered phase disappear, and 
anomalous exponents are confined to the immediate vicinity of the critical point. Note 
that according to (3.1a) m(£) diverges for e = 2, too, which implies that the above 
mentioned inconsistency when assuming a phase with spontaneously broken continuous 
symmetry for n > 2 at two dimensions, does not appear in the Ising case. Of course, 
contrary to the systems where the Mermin- Wagner- Ho henberg theorem [22,23] applies, 
the two-dimensional Ising model displays a second-order phase transition at a finite 
temperature. 

Fig. la shows the numerical flow diagram for the effective coupling u e ff(£) vs. the 
quantity m(£) 2 /[l + m(£) 2 ]; compare with Fig. 4 of Ref.[ll]. Again, for m(l) < 0.1 one 
finds scaling behavior with respect to the variable x of Eq.(2.23), which is not surprising 
as this is essentially a property of the critical point. Hence for temperatures not too 
far away from T c universal crossover features can be studied, see Fig. lb (and compare 
Fig. 5 of Ref.fll]), starting from the Ising fixed point 

«J = y (3.4a) 
C/ = -2 + e (3.46) 

and terminating at the asymptotically stable Gaussian fixed point with 

u* G = oo , u *ff = (3.5a) 
Cg = "2 • (3-56) 

We illustrate the ensuing crossover within the q- and w-dependence of the correla- 
tion functions for the case of the static susceptibility. Using the numerical solutions for 
u(£) and m(£), and inserting n = 1 into Eq.(5.11b) of Ref.[ll], one arrives at the one- 
loop cumulant (C). The corresponding vertex function (V) and the zero-loop expression 
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follow by expanding with respect to u or putting u = 0, respectively. For comparison 
of the different results the effective exponent 

2 " Veff = dlng (3.6) 

is plotted in Fig. 2. In all three cases shown, one finds similar features because grossly the 
flow of the parameter m(£) describes most of the crossover behavior. For large values of q 
one has ?y e // = (to one-loop order), before eventually the static susceptibility becomes 
finite for q^O. This is to be contrasted with the divergences induced by the transverse 
massless modes for the models with n > 2 (see Fig. 7 with y = 0). We finally remark 
that our interpretation for the minimum of the effective exponent 2 — Veff apparent for 
the time-dependent Heisenberg models at intermediate wavenumbers [11], namely that 
it is due to the freezing-out of the longitudinal fluctuations before the Goldstone modes 
become the predominant excitations, is quantitatively supported by Fig. 2. 



IV. THE INFLUENCE OF WEAK CUBIC ANISOTROPIES ON THE 
CRITICAL DYNAMICS OF MODELS A AND B BELOW T c 

A. Model and renormalization 

We now start our investigations of the influence of characteristic anisotropies with 
the discussion of the crossover properties when cubic anharmonicities have to be taken 
into account. Thus we add the term 

/n 
^5>o(^) 4 (4-1) 

' a=l 

to the Hamiltonian (2.1a). Its static properties above and below T c were examined 
by Wallace [18] and Aharony [24] on the basis of a 1/n-expansion, and by Ketley and 
Wallace [25] using the e-expansion. The phase diagram of the cubic model was described 
in more detail by Brezin, Le Guillou and Zinn-Justin [26], Nattermann and Trimper [27], 
Lyuksyutov and Pokrovsky [28] , Rudnick [29] , Iacobson and Amit [30] . For a summary 
of their main results we refer to Chap. 4 of Ref. [20] . We shall restrict ourselves to the case 
vq < 0, for then the spontaneous order parameter stays along one of the principal axes 
of the n-dimensional hypercube. For stability reasons we furthermore have to require 

b = u + v o > ; (4.2a) 
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for b = the transition becomes first-order [27-30]. Finally we are only interested in a 
small deviation from the (3(n)-symmetric case and hence tiny anisotropy parameters 



0<y = ^«l , 
fro 



(4.26) 



such that the continuous character of the transition is guaranteed [28] . 

Eq.(2.3) again constitutes the Langevin dynamics of the system, supplemented by 
the moments (2.4) for the stochastic forces. For the spontaneous order parameter of 
Eq.(2.9a), however, we write 

4 > o = \7-m Q (4.3) 
V b Q 

instead of (2.10), and the fluctuation-induced T c -shift as determined from Eq.(2.11) for 
mo = acquires an additional contribution oc v o 

2/e 



roc 



U + 2 A 1 

— — u Q A d + — v A d 
6e 2e 



(4.4) 



with respect to the isotropic case. 

Hence in the ordered phase one finds for the cubic Janssen-De Dominicis functional 
J = Jq + Jint + Jct, where the harmonic part explicitly reads 

M{^o}^o,Uo}^o]= / / $^A g a ^o(£^)^o(-<f, -^) + X q a a (q,uj)ao(-q,-uj) 

- 2^o(£ w ) iuJ + Xo q a (ml + q 2 ) 7Tq (-g, -a>) 



cr (q,u) iu + \ q a (ml + q 2 ) cr (-q,-<jj) 



while the interactions are given by 



(4.5a) 



OJi x 



X 



u o ^i) TToifo, uj 2 ) TTo (93, w 3 ) 7rJ(g 4 , w 4 ) 

a/3 

+ ^2 V ^0 (si > ^l) < (92, W 2 ) 7To (?3, W 3 ) 7Tq (g 4 , W 4 ) 
a 

+ WO 7fo (91,^0 ^O (92, W 2 ) 0-0(^3,^3) 0-0(^4,^4) 

a 

+ wo cro(9i,wi) ttq (9*2,^2) 7Tq (93,^3) cr (<f 4 ,u; 4 ) 
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+ fro croCsi, ^i) cr (<f 2 , w 2 ) cr (<f 3 , W3) ^o(?4, w 4 ) 



(4.56) 



l x °\ h~ m ° I I q " 6 ( s * ) 5 ( Wi ) x 

D V °0 Jqxqiqz \ , / \ i J 

^2 2 W TTo (?i > ^l) ^0 (?2, W 2 ) CTo(?3, W3) 

a 

+ J^tiO cro(5l,u>i) 7Tq (52,^2) TTo (93,w 3 ) + 3 6 cro(?l,Wi) 0-0(92,^2) <T (93,^3) 

a 

and the counter-term stemming from the equation of state (2.11) finally is 
JcT[{TT%},ao,{Tr%},a ] = -\ A / q a itg(q, u)ir%(-q, -u) + a (q, co)a (-q, -u)] 

J q J us „. 



A J-m o iff O (0,0) o a ,o 



(4.5c) 



The most important difference to the isotropic model, which of course emerges as the 
special case with vq = and 60 = uq, is the appearance of the transverse mass (remember 
v < 0) 

— 2 ~ v 2 1 A a\ 

m ° = 26o~ m ° ^ ' 

in the harmonic functional (4.5a). The basic elements of the perturbation theory to be 

constructed from Eq.(4.5) and their graphic representations are depicted in Fig. 3. 

The ultraviolet divergences of the dynamic correlation functions are multiplicatively 

renormalized by introducing the Z-factors 



2 ry-\ 2 -2 
2 ry- 1 2 -2 

m = z w m o V 
b = Z' 1 b A d fj>- 

V = Z' 1 Vq A d \l 



(4.7a) 
(4.76) 
(4.7c) 
(4.7d) 



in analogy to (2.15). In the framework of a one-loop theory we may omit the field renor- 
malizations and hence the Z-factor for the Onsager coefficient A. The renormalization 
constants introduced in (4.7) are not independent; e.g. from the definition (4.3) one 
easily derives 

Z b = Z m Z a , (4.8a) 
with Z a = 1 to one-loop order. Similarly, (4.6) leads to 



Zb Zjn — Z v Z m 



(4.86) 
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The renormalization constants (4.7a,b) are conveniently determined by rendering 
the quantities (c^)"^ 2 F^ a (q, u) and (c^)^ 2 r^ 7r (g, uj), respectively, finite at the nor- 
malization point q = fx, ui = 0. Using the explicit expressions for the two-point vertex 
functions listed in the Appendix, the one-loop results are 



Z m = Z b = 1 + 



n-1 (6 - v ) : 



6e 



bo 



3_ b A d fi~ 



e/2 



(4.9a) 



1 + 



1 4b -v v Q A d /j~ 



+ 



1 46 - vq b A d [i 



2e2b + v Q (i + ml/^y/ 2 e2bo + v (i + m 2 /n 2 ) e/2 

Obviously, for v = (4.9a) reduces to Eq.(4.11c) of Ref.[ll] for the isotropic model 
The C- and /3-functions following from (4.9) read 



(4.96) 



Cr, 



d_ 

d 



In 



m 



-2 + 



n — 1 



(b-v) 2 /b 



+ 



o mt 



b = b 



6 (1- frm 2 ) l +t/ 2 2(l + m 2 ) 1 + e / 2 



-e + 



n — 1 



(6- v) 2 /6 



+ 



A, = A» 



d_ 



v = v 



-e + 



6 (l-^m 2 )^/ 2 2(l + m 2 ) 1 + e / 2 

Ab-v ( v/2 b 

+ 



(4.10a) 

(4.106) 
(4.10c) 



2b + v ^m 2 )^/ 2 (l + m 2 )^/ 2 

and the flow-dependent couplings are given by the following set of ordinary diffential 
equations 



£ 



dm(£) 
d£ 



1 



m(e)Cm(£) 



dv(£) 



£ 



d£ 



(4.11a) 
(4.116) 
(4.11c) 



with the usual initial conditions m(l) = m, 6(1) = 6 and v(l) = v. The renormalization 
group equations for the two-point cumulants are finally solved by [compare Eq.(2.18b)] 

G% w (n, A, m, 6, v, {q}, M) = G%M A, m(£), b(£), v(£), {q/fi£}, {u/^£) 2 }) . (4.12) 



B. Discussion of the flow equations and fixed points 

The common zeros of the (3- functions (4.10b,c) yield the fixed points governing the 
universal critical behavior 



P b (b*,v*,m) = = p v (b*,v*,m) 



(4.13) 
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Their regions of stability can be inferred by investigation of the matrix 

O = 



d(3 h /db dPb/dv 
d(3 v /db d(3 v /dv 



(4.14) 

b=b* , v—v* 



(with fixed m): if both eigenvalues of are positive, the fixed point is (infrared-) stable, 
one positive and one negative eigenvalue correspond to a saddlepoint, and two negative 
eigenvalues lead to instability. 

We begin with the case m = 0, i.e.: the properties of the critical point T = T c . 
Besides the Gaussian fixed point b* G = v G = 0, which is unstable below d c = 4 because 
both eigenvalues of are equal to — e, one finds three non-trivial solutions of Eq.(4.13), 
namely (i) the isotropic Heisenberg fixed point (see Eq. (2.20a) of Section II. B) 

b * H = ^fs ' V * H = ° (4 ' 15a) 

/l 2(n-l) \ 

Vh= ( )e , (4.156) 

^ U n+8 / 

(ii) the Ising fixed point [compare Eq.(3.4a)] 

b*j = v*j = y (4.16a) 

^/=(^_ 1 4 _°i) e ( 4 - 166 ) 

with saddlepoint-type behavior, and finally (iii) the cubic fixed point 

2e n - 1 _ 2en-4 
^4 = t; , ^4 = -= (4.17a) 

A 3 n A 3 n v ' 

n ( 5n-8 2(n-l)\ e 

where has the eigenvalues e and (n — 4) e/3n [20 (and references therein)]. Therefore 
the situations n > n c and n < n c have to be distinguished, where to one-loop order 
n c = 4 (in more elaborate calculations one finds that n c is between 3 and 4) [25-27]. 
For n > n c the cubic fixed point (4.17a) is stable, and hence the anisotropy (4.1) 
constitutes a relevant perturbation of the 0(n)-symmetric Hamiltonian (2.1); to the 
contrary, in the case n < n c the Heisenberg fixed point (4.15) becomes asymptotically 
stable, which means that the fluctuations dynamically restore the rotational symmetry. 
(An interpretation for this remarkable difference might be that in high-dimensional 
spaces most of the "volume" is concentrated in the vicinity of the surface, while for 
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low dimensionality n < n c the detailed shape of the surface in configuration space is 
less relevant and the "edges" are effectively smoothed by the fluctuations.) For n = 4 
Eqs.(4.15) and (4.17) are identical, and of course for n = 1 only the Ising fixed point 
(4.16) remains and becomes stable, compare Section III. In the limit n — > oo (4.16a) 
and (4.17a) coincide. 

From (4.15)-(4.17) we can deduce the qualitative features of the flow diagrams, 
which are displayed in Fig. 4 for the typical cases n = 2 and n = 8, respectively. In 
Fig. 4a the attractivity of the Heisenberg fixed point "H" (4.15) is evident, whereas in 
Fig. 4b the anisotropic cubic fixed point "A" (4.17) dominates the flow of the couplings. 
Those trajectories which are not located in the regime of the respective stable fixed point 
eventually cross the stability boundary 6 = and the corresponding physical system 
will then display a first-order phase transition [28-30,20]. 

We can now turn to the phase with spontaneously broken (discrete cubic) symmetry, 
i.e.: the case m — > oo. Apart from the (for e > 0) unstable Gaussian fixed point 
b* G = Vq = there is now the isotropic coexistence fixed point [see (2.21a)] 

b*c=^T[ ' v c = ° ( 4 - 18 «) 

n C =(l \)e , (4.186) 



x _1 , 

which, however, in contrast to the isotropic model behaves like a saddlepoint in param- 
eter space. As long as v ^ 0, we have asymptotically m(£) — > oo. Hence from (4.10) 
and (4.11) we conclude 

m{£) oc r 1 (4.19a) 

b(£) oc £~ e (4.196) 

v(£) oc £~ e , (4.19c) 

i.e.: the mass parameter and anharmonic couplings diverge in a manner completely 
analogous to what we have previously discussed for the case n = 1. Of course, a 
crossover to an asymptotically Gaussian theory should in fact be expected because the 
fluctuations are bound to cease eventually. This can be confirmed by studying the 
effective couplings 

KffV) = M , JU» « f (4.20a) 
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compare Section III and Fig. 5. In particular, we want to point out that in the presence 
of cubic anharmonic terms the Mermin-Wagner-Hohenberg theorem [22,23] does not 
apply and hence long-range order is possible in two dimensions, which is reflected in 
the ^-independent divergence of the mass parameter according to (4.19a) instead of the 
anomalous dimension (2.21b) of the isotropic model. 

The fact that (to one- loop order) for n < 4 the 0(n)-symmetry is dynamically 
restored at the critical point turns out to be most important for the flow of the effective 
coupling b e ff(£) and hence the static and dynamic susceptibilities, at least for low 
initial values of the anisotropy parameter y of Eq.(4.2b). Ideally the system would 
acquire v* H = at T c leading to pure coexistence behavior as described in Section 
II. However, it is unstable towards any small disturbances in the indirection followed 
by a crossover into the asymptotically Gaussian model. We may study this situation 
quantitatively by changing the initial value y(l) in the numerical solution of the flow 
equations (4.10), (4.11). In Fig. 5 the flow of the effective coupling (4.20a) vs. the scaling 
variable (2.23) is depicted for several values of y(l); for comparison the crossover from 
the Heisenberg to Lawrie's coexistence fixed point in the (9(n)-symmetric model is 
plotted, too. For very low anisotropy b e ff(£) rises towards the coexistence fixed point 
(4.18), but finally approaches the Gaussian limit [see Eq.(3.5)]. The smaller y(l), the 
higher is the maximum of b e f / (£) , and the later the crossover to b* e ^ — takes place. In 
fact the approach to b e f /(£)—> could be universally rescaled by using the " Gaussian 
scaling variable" £/m(l) instead of (2.23). 

The large attracting region of the isotropic coexistence fixed point (4.18) becomes 
obvious from the qualitative flow diagram in Fig. 6a (with n = 2; compare Fig. 4a). For 
n > 4, on the other hand, for the allowed values v < the Gaussian limit is immediately 
assumed, without any relics of the coexistence region, see Fig. 6b. Even for very low 
values of y(l) no quasi- Goldstone modes are of importance. We finally remark that the 
above considerations, especially for the critical region m = 0, are inferred from the one- 
loop approximation only. Thus by taking into account higher orders of the perturbation 
expansion modifications might be necessary. In particular, the number pf components 
n c below which the Heisenberg fixed point is the attractive critical fixed point is shifted. 
It should be clear, though, how the properties of the critical point imply the asymptotic 
features for T < T c and q — > 0, uj — > 0. 
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C. Dynamic response and correlation functions 

With the help of Eq.(4.12) we can now evaluate the dynamic two-point cumulants 
to one-loop order from the diagrams listed in the Appendix. We furthermore apply 
the general matching condition (2.24) designed to connect the low-wavenumber and 
low-frequency behavior with the asymptotic limit I — > where the stable fixed points 
are approached. In the presence of cubic anisotropies the perturbation series cannot 
be truncated at the one-loop level, however, as opposed to the isotropic case. For this 
reason, and to avoid somewhat even lengthier expressions, we quote the results for 
the one-loop vertex functions which may be easily transformed into the cumulants for 
comparison with the 0(n)-symmetric models [11]. Specializing (2.24) to u> = 0, using 
Eq.(4.9) and the abbreviation 



W = ,, 1 + M^ + M^ , ( , 21) 



the inverse static susceptibilities at d = 3 (e = 1) read (see Eq.(5.5) of Ref.fll]) 



_ [b(£)-v(£)] 2 2 1- 2gg ™W 1 + m m(£) 

6b(£) m[i) r CSm S(£) +arCSm S(£) 

and 



1 + m(£) 2 Z m (£) - - b(£) m(£) 2 arcsin - 1 (4.226) 
2 i/l + 4m(l) 2 



n - 1 [fe(£) - v(£)] 2 fe . 2 . 1 
J m(ty arcsin 



6 6(£) w " /-, | -2«(£) 



£=q/ fJ , 



Again, their features are most conveniently explored by introducing effective expo- 
nents 



»/-_ . (4-23) 



^T/L _ gln X T/L (9,0) 
91ng 

which are displayed in Fig. 7 for n = 2 and the same values for the anisotropy parameter 
y(l) that have been used in Fig. 5, including y(l) = 0. (In order to allow for a most direct 
comparison with the isotropic model, here the two-point cumulants are plotted.) The 
transverse exponent (Fig. 7a) shows how the quasi- Goldstone modes for y(l) ^ finally 
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acquire a mass and hence the transverse susceptibility becomes finite for q — > 0. If we 
had plotted 2 — ryj^ against q/fj,m(l) instead of the critical scaling variable following 
from (2.23), we would have found a universal crossover signature very similar to the one 
depicted in Fig. 2 for the one-component case. For very low anisotropies, on the other 
hand, coexistence-like anomalies are clearly noticeable in the wavevector-dependence 
of the longitudinal response function within an intermediate (/-region the size of which 
depends on y(l), if n < n c . For y(l) < 10~ 4 , the pronounced minimum of the effective 
exponent 2 — is then still apparent, while for larger anisotropy the flow becomes 
much more similar to Fig. 2 for the direct crossover to the Gaussian model. Now turning 
to n > n c , according to the preceding Subsection, no remnants of the Goldstone modes 
are to be expected. 

As in the isotropic case, the dynamic response functions may be explicitly written 
down for arbitrary dimensions 2 < d < 4 in the limit q — > 0. For model A [1] with 
non-conserved order parameter one finds (compare Eq.(5.11) of Ref.fll]) 



2b(£) + v(£) 



-v{£) 



m(£) 2 Z w (£) 



x 



+ 



+ 



6e 
m(£) 2 



1 - 



% — 



2b{£)+v{£) 
2b(£) 

m(£) 2 



m {ey 



2b(£) 

—v(i) x 
2b{£) 

-v(£) 



2-e m-v(£)] 2 



2b(£) 



m(£f 



i + 



2b(£)+v(e) 
2b(£) 



m(£y 



m ® - 3eb(£) X 

1_e/2 f-i , 2b( £)-v(£) 
4b{£) 

l-e/2' 



m(ey 



I -i 2b{£)-v{£) . n . 2 



(4.24a) 



l-e/2' 



£=^/u/\^ 



Xl l (0, to) = y 2 £ 2 i-i + m{l) Z m (£) + ^ b(£) m{£f 
(n-l)i [b(£)-v(^ 2 



m{£) 



2-e 



+ m(ey 



+ 



X 



3e 



b(£) 



m(£Yx 



l-e/2' 



(4.246) 



2b(£) 



m{£)' 



l-e/2 



—i -v(£) ...o 



l-e/2' 



£=^/co/\^ 



The logarithmic derivatives of the inverse real part of the longitudinal susceptibility 
and of the inverse longitudinal correlation function obtained with the help of (2.8) are 
depicted for n = 2 and e = 1 in Figs. 8a and 8b, respectively. For the relevant situation 



22 



n < 4, their general features are similar to that of the static effective exponent 2 — Veff'i 
however, the crossover to the asymptotic Gaussian behavior takes place for values of 
the dynamic scaling variable about an order of magnitude lower than in the static case, 
and the influence of the cubic anisotropies is considerably reduced. 

If the order parameter field is conserved (model B according to Hohenberg and 
Halperin [1]), due to the diffusion pole the ratio u/q 2 has to be kept fixed when the 
limit q — > is performed. Then merely the zero-loop contributions (to be supplied with 
the full one- loop flow equations) persist [11]: 



Xl 1 (0,u/q 2 ) = fie 2 \-i + m(£) 2 



m(ff 



(4.25a) 
(4.256) 



e=y/w/\q* 

The ensuing crossover to the Gaussian theory with pronounced remnants of coexistence 
anomalies again is shown in Fig. 9, for the case n = 2 and d = 3. 

For a comparison with the results (5.13) of Ref.fll] we finally quote the complete 
q- and w-dependent response functions for model A at three dimensions 

2 6(f) +v(£) 



-1 / -> \ 2 IU 2 ifl 



2b(£) 



m(£) 2 Z w (£) + 



6 



1 



[b(£)-v(£)} 2 m{£f 



+ arcsin 



6 6(f) q/n 
iuj/X 



arcsin 



q 2 -iu/X 2b(£)+v(£) 



fi 2 £ 2 S+(£) 2b(£)S+(£) 



m{£)' 



/j 2 £ 2 S+(£) 2b(£)S+(£) 
+ arcsin 



+ arcsin 



2 6(f) 



2b(£) + v(£) ,,, 2 



+ 



H 2 £ 2 S_(£) 
2b(£) + v(£) 



q 2 — iu/X 

ix 2 £ 2 S-(£) ' 2b(£)S-(£) 



2b(£)S_(£) 
m{£) 2 




and 



X - L \q^)=q 2 -^+^ 2 £ 2 



m{£) 2 Z m {£) 



n - 1 [6(f) - v 



m(£) 2 ( 



6 6(f) q/ftf V 

^2 



arcsm 



iu/X 



_ 3 b{£)m{£) 2 ( 
~ 2 q//i£ 



arcsin 



H 2 £ 2 S ± (£) 
iu/X 



+ arcsin 



V 2 £ 2 S l{ (£) 



+ arcsm 



q 2 — iu/X 
lx 2 £ 2 S L {£) 

iu/X 



(4.266) 



V 2 £ 2 Sn(£) 



Here the abbreviations 

S±(£) = 



% u/X ± 2b(t) + v(£) 



jx 2 £ 2 



26(f) 



+ 



g 2 ~ 2 *^M + nz)-v(e) m{£)2 



H 2 £ 2 



6(f) 



H 2 £ 2 
(4.27a) 
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and 




^xY + (l^ + = im m(e) ,)jf (4 , 76) 



n 2 i 2 ) \ ^ 2 i 2 b(£) v ' ) ^i 2 

iuj/\\ 2 {q 2 — 2iuo/\ Q 2 ,. . -> 

1 +[—. T7T-+4mW J 4u (4.27c) 



fj?e 2 J \ v 2 i 2 K ' ) [i 2 t 2 

were used, and the matching condition 



2 



g 2 z u; 



(4.28) 



// 2 A /U 2 

has to be inserted. Investigation of the poles of (4.26) shows that the damping of both 
the transverse and longitudinal modes asymptotically is independent of q. Similarly, in 
the case of diffusive relaxation (model B) the dispersion relation becomes % uj oc q 2 . The 
anomalous line shape (2.27) of the longitudinal correlation function still appears in an 
intermediate range, if n < 4 and y(l) -C 1; finally it develops into the regular Lorentz 
form. 

The graphs discussed above refer to the situation where the number of components 
is low, n < n c . For larger n, the cubic anisotropy immediately destroys the massless 
character of the transverse excitations and hence no coexistence- type anomalies persist. 
The theory is generally restricted to dimensionality 2 < d < 4. At the upper critical 
dimension d c = 4, instead of the power laws logarithmic corrections come into play. 
Although excluded from the scope of the formalism presented here, in the limit d — > 2 
reasonable statements may be extracted from the general features of the flow equations, 
in the sense that the results are compatible with exact theorems on the existence of 
long-range order [22,23]. In Table II we summarize the qualitative behavior under the 
influence of small cubic anisotropies. 



V. THE INFLUENCE OF THE DIPOLAR INTERACTION ON THE 
DYNAMICS OF MODEL A IN THE ORDERED PHASE 

A. Dipolar propagators in the ordered phase and 
crossover behavior at the critical point 

The second very interesting generic anisotropy we want to study is the dipol-dipol 
interaction. Although the phase transitions, e.g. in magnetic systems, are driven by 
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the short-range exchange interactions, long-range dipolar forces become increasingly im- 
portant on approaching the critical point where the correlation length diverges. Within 
the static renormalization group theory, the influence of dipolar terms has been exam- 
ined by Fisher and Aharony [31], and was systematically treated together with other 
kinds of anisotropies by Aharony [24]. The interplay between the isotropic Heisenberg 
fixed point and an additional asympotically stable dipolar fixed point produces remark- 
able crossover phenomena, which have been intensively investigated in the literature; 
we quote here the works of Bruce and Aharony [32], Nattermann and Trimper [33], 
Bruce, Kosterlitz and Nelson [34] , and later of Santos [35] , Kogon and Bruce [36] , which 
are all based on the renormalization group theory, and expecially the recent papers by 
Frey and Schwabl [37] who applied the extended renormalization scheme of Amit and 
Goldschmidt [14] to the dipolar crossover problem. 

Taking into account only those terms which are relevant with respect to the renor- 
malization group [31], the starting point of our considerations is the following supple- 
ment to the static Ginzburg-Landau functional (2.1) 



Hence the dipol-dipol interaction hence leads to a term proportional to the coupling 
parameter go which is anisotropic in momentum space. We shall treat the general 
situation n ^ d, although in part of our study we shall restrict ourselves to the widely 
considered case of equal dimensions in momentum and order parameter space. With 
the aid of the projectors [37] 




min(d,n) 



(5.1) 



<*,P=1 




S a ~ 



9a q& 



(5.2a) 




(5.26) 



fulfilling the relations 




(5.3a) 



d 




(5.36) 



7=1 



d 




(5.3c) 



7=1 
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we may decompose the Hamiltonian into different contributions corresponding to the 
transverse (T) and longitudinal parts (L) with respect to the wavevector q 

min(d,n) 



#o[{0o a }]=/ o E Wo + q 2 )P^ + (r + go + q 2 )P^Uo(Q)^o(-Q) 

J *V a,/3=l 



(5.4) 



E toMM-fi 

a=min(d,n) + l 



From (5.4) it is apparent that the fluctuations in the longitudinal sector acquire an 
additional mass term oc go due to the dipolar forces. 

The static critical properties of any dipolar magnet are described by the Hamilto- 
nian (5.1) irrespective of the detailed dynamics. For the latter we remark that in the 
presence of dipolar forces the magnetization is no more a conserved quantity. Hence 
the relaxation process would be that of model A [1], Eq.(2.3) with a = 0. However, 
mode-coupling terms are most important for the dynamic properties of ferromagnets and 
antiferromagnets. Yet for the sake of simplicity, and because the focus of our intention 
is rather the influence of the anisotropies on the fluctuations, we shall drop the mode- 
coupling vertex and restrict our investigations to the time-dependent Ginzburg-Landau 
model A. 

In order to avoid unnecessary complications, we start with the important special 
situation where 

n = d = 4-e , (5.5) 

and shall come back to the general case later (Subsection C). Below the transition 
temperature, we may still use the parametrization (2.9), (2.10) for the spontaneous order 
parameter. The Janssen-De Dominies functional (2.5b) can then be split according to 
J = Jq + Jq + Ji nt + Jar, with the harmonic contributions 

^o T [R Q }^o,{<},^o] = (5.6a) 
= / /[E (^*o(?>")^(-£-")-*o(£w) [^ + Aog 2 ] 4(-q,-uj))p^ 

+ ^\ ao(q,u)a (-q, -u) - a (q,uj) iu + X (ml +Q 2 )] <r (-q, -uj)^j 

and 



Jo[{ 7t o}^ a o,{ n o}^ a o] = 



(5.66) 
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L 

a(3 



L 

nn 



= Yl ( A °^o -n%(q,uj) i u + X (g + q 2 ) ir${-q, -to) J P t 

+ (a o- (q,u) 0- o (-q,-u) - a (q,cu) iu + X (ml + g + q 2 ) a (-q, P % 

a 

denoting the transverse and longitudinal parts with respect to the wavevector q. The 

interaction term Ji nt is identical to that of the isotropic model A and follows from (4.5b) 

for vq = and bo = uq. Similarly, the counterterm Jqt looks like (4.5c), with A replaced 

by the analogous quantity A. In the transverse part (5.6a) of the dynamic functional the 

7r-modes are massless and in harmonic approximation independent of the cx-fluctuation. 

The longitudinal 7r-modes, on the other hand, have lost their Goldstone character due 

to the dipolar interaction, and furthermore couple to the cr-field (oc Ao<7o)- We also 

remark that the operators P^p and P^ are no proper projectors in the subspace of the 

n — 1 fluctuations transverse with respect to the direction of the order parameter. 

The evaluation of the harmonic propagators is thus more difficult than in the 

isotropic case. It is very convenient to collect the response and correlation fields in 

a 2 n-dimensional vector according to ifo = (7Tq, . . . , 7Tq _1 , ao, 7Tq, . . . , tyq -1 , ao). Then 

Jo = Jq + Jq can be written as 

1 



ip (q,uj) A(q,oj)ip (-q,-uj) 



qJ lo 



where the harmonic coupling matrix separates into n x n-submatrices 

D G 



Here D is diagonal 



D = 








-2A 



and the symmetric matrix G has the form 

G = 

with 



B v 
v T 6 



B a( 3(q, u) = (i lo + A q 2 ) 5 a p + A go 



v a (q) = X g 



q a q r 



q 2 

b(q, to) = i uj + A (m + q 2 ) + X Q g 



(5.7a) 
(5.76) 

(5.7c) 

(5.7d) 

(5.8a) 
(5.86) 

(5.8c) 
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The (n — 1) x (n — l)-matrix B is easily inverted: 



1 



The inverse of (? then follows as 



AoS'o 



«w + A ((/o + q 2 ) - ^ogoql/q 2 q 2 



(5.9) 



G -1 = 



with 



C u 

T 

u c 



1 



(5.10a) 



b-v T B~ 1 v 
u = —cB~ x v 

C = B' 1 (l-vu T ) 



Finally one finds 



G* _1 

G-i -G _1 DG* _1 



(5.106) 
(5.10c) 
(5.10d) 

(5.11) 



from which the propagators of the theory may be obtained (see Ref.[ll,12]). The zero- 
loop response propagators read 



Go^a a (q, u>) = 
G oa<7 (q,<jj) = 
here the abbreviation 
s(q,u) -- 



1 



-i iv + Xo q 2 
-1 



c q^qp x 

o a (3 n-s(q,uj) 

q a q n 



+ X (ml + q 2 ) q 2 



1 



-iu + A (m5 + q 2 ) 



1 _<ti -iuj + X q 2 



q 2 -iu> + Xo(ml + q 2 ) 



s(q,u) 



(5.12a) 
(5.126) 
; (5.12c) 



Ao go 



-ioj + Xo (g + q 2 ) - A go _ lLJ+x ° o( ^ +q2) q 2 Jq 2 
has been introduced. Note that in the static limit uj — > they acquire the form 



(5.12a 1 ) 



n f-\ 1 (x MP 1 

Go a 0{q) = ~~z : — o d «/3 ~~ 90 



r a _|_ ^,2 



"o + 9 2 ?2 + #° 



where 



9a 
r a _|_ ? 2 



(5.13a) 



(5.136) 
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which is precisely the result quite generally derived by Aharony [24] for anisotropic 
exchange accompanied by dipolar interactions, if Tq = for a = l,...,n — 1 and 
Tq = ttiq is inserted. 

The rather complex structure of the result (5.12) for the propagators of the har- 
monic theory renders a complete treatment of the crossover regime a very streneous 
task. Hence we shall restrict our discussion to the coexistence limit mo — > oo only, 
which suffices for the questions we have in mind. Before embarking on that, however, 
we are going to discuss very briefly the situation at the critical point. 

With m = 0, the response and correlation propagators at T c read 

=_ iu l Xoq2 P lp + — + xligo + a 2 ) P ^ 



2 A 2 A 

G ^ (q, U ) = ^ - ^ q4 Pip + w2 + A 2 J q + q2) 2 P tp ■ ( 5 - 14& ) 



For T 7^ T c , one may reintroduce the temperature dependence ro (see Ref.[37]). The 
T c -shift can be implicitly determined to first order in uq, e.g. by using the equation of 
state (2.11), as the solution of the equation 

2/(2-6) 

' O /y-t Q 1 

."0.1/ 



9o = r 0c 



rj 2n 3) - ! 



(5.15) 



Compared to the isotropic case, the transition temperature is increased due to the 
long-range order dipol-dipol interaction which favors condensation [31]. 

For the renormalization in the critical region conceptually different Z-factors for 
the transverse and longitudinal field components, respectively, have to be defined [37], 
namely 

r = zf £ p£ % + zf Yl p «p « ( 5 - 16 «) 

& p 

p p 

The renormalized counterparts for the timescale Ao and the coupling constants go 
are conventionally introduced as 

\ = Z- 1 X (5.17a) 
u = Z~ x u A d /i~ e (5.176) 
g = Z; 1 g ^- 2 . (5.17c) 
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As a consequence of the fluctuation-dissipation theorem (2.8) and the partition of the 
propagators into transverse and longitudinal components according to (5.14) the follow- 
ing relation between Z\ and the field renormalizations can be derived 

Zrp -\~ Z . . 

A= (z f Zt)-v 2 + (z l z l )-w ' (5 ' 18) 

which for Zj, = Z^ and Z? = Zl (valid for vanishing go) reduces to the well-known 
expression for the isotropic relaxational models [12,11]. 

In the vicinity of T c , interesting crossover phenomena occur which are induced by 
the second lengthscale 1 / ^/g~o originating from the dipolar term. An elegant description 
of these phenomena may be obtained by applying Amit and Goldschmidt's general 
method [14]; hereby it appears that there are no ultraviolet divergences oc go, and hence 
one concludes [37] 

Z; X = Z L . (5.19) 
To one-loop order the field renormalizations vanish 

Z f = Z L = Z T = Z L = l , (5.20a) 

and therefore 

Z x = 1 (5.206) 

Z g = 1 . (5.20c) 

For the single non-trivial Z-factor one finds by investigation of the four-point vertex 
function [37] 

n 2 + 18n + 116 _ e n 2 -6n-76 u A d ^ 

Z u = 1 ^ 771 u A d ii — — (5.20d) 

144e 144e (1 + g Q /^t 2 ) e / 2 

(here we have slightly modified the results from Ref. [37] by not explicitly inserting n = 4 
and refraining from the e-expansion) . 

Taking advantage of the renormalization group equation we introduce flow-depen- 
dent couplings according to 

Z dj ^=9WtS) (5.21a) 
t^p-=P*W , (5-216) 
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with the usual initial conditions g(l) = g and u(l) = u, where 



Cg = V 



d_ 



In 



-2 



o 9o 



and 



R 8 
Pu= 



u = u 



-e + 



n 2 + 18n + 116 



n 2 — 6 n — 76 



u 



u 



(5.22a) 



(5.226) 



144 144 (l + g) 1 +*/ 2 _ 

are Wilson's flow functions to one-loop order. For vanishing dipolar interaction g = 
one rediscovers the isotropic Heisenberg fixed point (2.20a), as the zero of (5.22b). For 
any finite value of #(1), however, the dipolar strength diverges as 



g(£) = gt 



(5.23) 



which ultimately justifies the corresponding term in (5.20d) and (5.22b). Near the criti- 
cal point where the correlation length becomes macroscopic, the dipol-dipol interaction 
hence is a relevant perturbation. The asymptotic behavior is then no more governed by 
the Heisenberg fixed point but by the dipolar fixed point [31] 

144 e 



U n = 



(5.24) 



D n 2 + 18n + 116 

The ensuing crossover features have been intensively studied by several authors [32-37] . 
For the possible values n = 2, 3, 4 one finds u* D > u* H , which accounts for the fact that 
a certain amount of the fluctuations is effectively suppressed by the dipolar forces. 

B. Renormalization of the theory in the coexistence limit for n = d > 3 

We shall now come to the properties within the ordered phase, where, as we shall 
see shortly, the cases n = 2 and n = 3, 4 have to be distinguished. For n > 3 the 
asymptotic theory can be discussed in a manner completely analogous to Section II. B. 
Substituting the transformed longitudinal fields (2.13) (with a = 0) into (5.6) and the 
interaction terms, and performing the limit mo — > oo one finds for the Janssen-De 
Dominicis functional 



^oo[W } 5 °"0, {tTq }, Co] 



(5.25) 



/ / ( A °^o (<?»^o -^o(<?» iuJ + Xoq 2 \n%(-q,-uj) 



af3 
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i.e.: remarkably a harmonic theory again, which can be treated exactly. 

One immediate consequence of (5.25) are the relations r^°- Q -^ (q, uj) = 2X 5 a p 
and r^. Q7r(3 (cf, u) = B a p [see Eq.(5.8a)]. Hence the field renormalizations for the tt- 
fluctuations vanish; a similar reasoning applies for the longitudinal fields, and from 
(5.18) and (5.19) one concludes that Eqs.(5.20) are valid as exact results in the coex- 
istence limit. In particular, for i — > the parameter g(£) diverges as in (5.23), and we 
simply have to insert s^q, u>) = q 2 /(q 2 — q n ) into the 7r-propagators: 

1 



P T 



n a Tcf 3 



-iu + X q 



2 A pT 

2 _i_ \2 „4 «/3 



UJ 2 + Agg 

Note that of the two proper projectors [compare Eq.(5.3)] 

Q a Q(3 



P L - 



q 2 -ql 



q 2 -q 2 



(5.26a) 
(5.266) 

(5.27a) 
(5.276) 



only the transversal operator appears in (5.26). (In the case q = (0, . . . , 0, q n ) one 
has P af3 = S a f3 and P af3 = 0.) The fluctuations parallel to the wavevector q have 
disappeared, as well as the non-diagonal propagator (5.12b). 

As for the isotropic model, using (2.13), (5.25) and (5.26) one finds the following 
asymptotically exact results for the cr-correlation functions 



(<?» = 



1 



A mg 



1 + 



U 



1 



1 



2 Y^ + ^+k 2 



x 



(5.28a) 



x 



n-2- 



(q 2 - Ql) (k 2 - kl) - [qk - q n k r 



2 1 



^Re 



(q 2 - Ql) (k 2 -kl) - (qk- q n k r 



(5.286) 



2 2 

^ + k 2 -k 2 n ) - (qk-q n k^' 



Again, these expressions can be identified with the explicit results from first-order per- 
turbation theory for the two-point cumulants in the coexistence limit. As for g Q = 0, the 
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corresponding two-point vertex functions may be represented by a geometric series of 
the 7r-loops, see Fig. 2 of Ref.[ll]. The essential differences between (5.28) and Eq.(3.18) 
of Ref.fll] are (i) the dependence on the angle between spontaneous order parameter 
and wavevector as apparent from the second term in the integrands of (5.28), and (ii) 
the factor n — 2 within the isotropic term instead of n — 1. 

For the determination of the only non-trivial renormalization constant = 
one has to consider rg°- cr (0, 0). Hence Z m is independent of the anisotropies in (5.28), 
and one finds explicitly that as a consequence of the dipolar interactions the relevant 
effective number of massless modes is reduced from n — 1 to n — 2, because 



n — 2 

Z™ = Z™ = l + —u A d ^ 



The remaining non-zero flow functions read 



C = -2 



n-2 
6 



(5.29) 



(5.30a) 
(5.306) 
(5.30c) 



and as a solution of (3^{u*) = we find the stable dipolar coexistence fixed point 

6e 



U CD — 



n-2 



(5.31a) 



replacing (2.21a); however, the anomalous dimension of the mass parameter is not 
altered 

CoD = -2 + e • (5.316) 

We remark that as a consequence of Z m = Z u this is a general result for every non-trivial 
fixed point value u* ^ 0. Therefore the power laws characteristic of the coexistence 
anomalies are not changed, either. For example, the dynamic susceptibility will have 
the asymptotic form 



Xl(<?» oc 



X/j, 2 



-e/2 



(5.32) 



leading to the anomalous lineshape of the longitudinal correlation function and the q- 
and ^-divergences summarized in the center column of Table I. But the amplitudes 
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of the scaling functions (5.28) are reduced with respect to the isotropic case, showing 
certain characteristic angle dependences. 

In two special cases, namely if the wavevector is either parallel or perpendicular 
to the spontaneous order parameter, the ^-dependent contribution in the bracket of 
Eqs.(5.28a,b) vanishes, and the expressions (5.28) may be directly compared with the 
isotropic correlation functions for a = 0. Due to the renormalization-group invariant 
(2.22), asymptotically m(£) 2 = m(l) 2 (u*/u(l)) £~ 2+e , and hence the amplitude of the 
longitudinal susceptibility is proportional to 1/u*. Inserting the fixed point values 
(5.31a) and (2.21a), respectively, we therefore we find the following exact amplitude 
ratio of the longitudinal response functions in the dipolar and isotropic case 

Xl(q,v) dipolar _ n — 2 ^ ^ 

X.l(Qi ^) isotropic re 1 

For n = 3 this universal amplitude ratio is 1/2, in accord with the results of Pokrovsky 
[38] , and Toh and Gehring [39] , obtained in the framework of a spin- wave theory. If the 
wavevector q is neither parallel nor perpendicular to the magnetization, the amplitude 
of (5.32) will be further reduced. 

For n = d = 4 there will be logarithmic corrections instead of power-law singular- 
ities. As is apparent from (5.29)-(5.33), the situation at n = d = 2 requires a separate 
discussion. Afterwards an immediate generalization to the general case n ^ d will be 
possible. 

C. The cases n = d = 2 and n ^ d 

The form of the dipolar coexistence fixed point (5.31a), in comparison with (2.21a), 
suggests that the effective number of critical fluctuations is reduced to n — 2. Hence 
for n = 2 there are no more massless modes left that may lead to infrared singularities. 
Indeed, the transverse projector (5.27a) disappears in this case, and inserting n = 2 into 
(5.30) one arrives at the following asymptotic behavior for the ^-dependent parameters 

g{£) oc r 2 (5.34a) 
m(£) oc i~ x (5.346) 
u(£) oc £~ e . (5.34c) 

According to the preceding investigations in Section III and Section IV. B we may readily 
identify this as the signature of a crossover to a Gaussian theory. The fluctuations die 
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out on leaving the critical region T ps T c , the a-mode because of the formation of long- 
range order, and the 7r-excitation due to the dipolar coupling. The above discussion 
of the asymptotic theory thus does not apply for the two-component model, and the 
effective exponents characterizing the static susceptibilities, for example, will rather 
look similar to Fig. 2. 

Remarkably, there appears no inconsistency for the assumption of long-range order 
in two dimensions, as opposed to the isotropic model. This, again, is no contradiction 
to the Mermin-Wagner-Hohenberg theorem [22,23] which is based on the assumption of 
solely short-range isotropic interactions. We would like to emphasize once again that 
more complicated situations, as inhomogeneous ordering or topological excitations, are 
not within the scope of the present formalism. 

We shall now drop the restriction (5.5) and return to the general situation where 
the number of components n and the spatial dimension d may be different. Let us begin 
with the case d > n; the upper limits of the sums in (5.1) and (5.4) are then n and 
the second term in (5.4) is absent. Hence the structure of the asymptotic Janssen-De 
Dominicis functional (5.25) is not altered. We only have to be careful when inverting 
matrix B of Eq.(5.8a), resulting in 

1 / r. Xogo 



iuj + \ q 2 \ ap iu + \ (g + q 2 ) - \ g (I ~ Q 2 /q 2 ) q 2 
instead of (5.9); here we have used the notation 



9a qp 



(5.35a) 



Q= (?l,---,9n-l) 



(5.356) 



In the coexistence limit the transverse propagators read 

l / - q a qp 



(<?>) = 



—% u) + A q 2 
2A 



5 a [3 - 



Q 2 



ui 2 + \ 2 q 4 V* P Q 2 
and for the longitudinal correlation functions one arrives at 

1 A Un f 1 1 



— 1 + ^ 
Xom 2 { 6 Jkfi_k) 2 ^ + ^ + k 2 



x 



X 



72-2 



Q 2 K 2 - (QK) 

(^ + k 2 ) 2 -(qk)\ 



(5.36a) 
(5.366) 



(5.37a) 
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Q 2 K 2 - (Q 



which obviously generalizes Eq.(5.28). 

In particular the renormalization constants are identical to the previous results; 
therefore for d = 4 and n = 3 there will be logarithmic corrections characteristic for the 
upper critical dimension, and for d = 3, 4 and n = 2 a crossover to a Gaussian behavior 
will be found. Thus for d > n there are no specifically interesting modifications with 
respect to the cases under discussion before. 

In the opposite situation, n > d, it is crucial whether the spontaneous order is 
created in one of the directions that are affected by the dipolar interaction, or not. One 
would expect the former situation to apply, because the dipolar forces favor conden- 
sation. Yet if there are further anisotropies present, even transitions between different 
orientations e.g. of the magnetization of two-dimensional Heisenberg ferromagnets may 
occur [40]; consequently we shall take both possibilities into account. The modified 
matrix keeps the form (5.35a), albeit with 

Q= (q n -d+i,---,q n -i) , (5.38a) 
if the a-modes are subject to the dipol-dipol forces, and 

Q=(qi,...,q d ) , (5.386) 

if they are not. In the first case there will be d — 1 transverse propagators of the type 
(5.36) with a, f3 = n — d + 1, . . . , n — 1, and n — d propagators 

Go°W (q, (J) = I <W (5.39a) 
— i to + Ao q 

G^*M") = 2^4 <W (5-3%) 

with a, [3 = 1, . . . , n — d that are equivalent to those of the isotropic model. If the a-field 
is not affected by the dipolar interaction, the first d transverse modes have the form 
(5.36), while the remaining n — d — 1 propagators are given by Eq.(5.39). 
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The cr-correlations and hence the longitudinal (with respect to the order param- 
eter) response functions look like (5.37), where one has to insert (5.38a) or (5.38b), 
respectively. The relevant asymptotic Z-factor is again (5.29), and for e > the dipolar 
coexistence fixed point determines the emerging coexistence anomalies. The leading 
term (in the limit q — > and u — > 0) of the longitudinal susceptibility is given by 
(5.32), accompanied by an angle-dependent scaling function which assumes a maximum 
if Q = 0, and is otherwise smaller than unity. 

Thus the results of Subsection B can be taken over quite easily to the general case 
with arbitrary number of components n > 2 and dimension 2 < d < 4. The qualitative 
features are summarized in Table III. 

VI. SUMMARY AND DISCUSSION 

For isotropic models, n — 1 massless Goldstone modes appear in the phase with 
spontaneously broken (3(n)-symmetry. Their origin is that no special direction of the 
spontaneous order parameter is preferred by the static Hamiltonian (2.1); hence an 
infinitesimal rotation of costs no (free) energy. These gapless excitations induce so- 
called coexistence anomalies in the entire ordered phase, see Table I [11]. Within our 
formalism [11], which is based on the crossover theory of Amit and Goldschmidt [14] and 
generalizes Lawrie's work [8] to dynamical problems, their analytical form can be traced 
back to the renormalization group fixed point Uq (2.21a), or rather the accompanying 
anomalous dimension (2.21b). Remarkably, the coexistence limit (T < T c , h = 0, 
q — > 0, u — > 0) can be exactly represented by a one-loop theory for the two-point 
cumulants or the leading order (spherical limit) of the 1/n-expansion for the two-point 
vertex functions. 

For the application to real systems, in particular to solids displaying structural or 
magnetic phase transitions, quite naturally the question arises whether certain anisotro- 
pic terms in the free-energy functional will completely destroy the coexistence anomalies, 
or if, and then under which circumstances, remnants of the Goldstone excitations will 
still affect the static and dynamic properties. In order to elucidate this issue, we have 
studied two characteristic cases, namely (i) cubic anharmonicity, and (ii) dipol-dipol 
interaction. Generically, three possible situations are conceivable: (a) the anisotropy 
immediately suppresses the transverse fluctuations, as will be the case for uniaxial dipo- 
lar forces or any mechanism directly introducing transverse mass terms, and also for 
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the cubic terms provided n > n c . However, (b) there may also be a quasi-Goldstone 
behavior in an intermediate range within the crossover from the critical to the Gaus- 
sian regime, especially if near T c fluctuations tend to restore the (9(n)-symmetry, an 
example of which is the cubic model with n < n c [25-27]. Even more subtle is the 
dipolar interaction, where (c) although the model is no more invariant with respect to 
a continuous symmetry transformation, not all the transverse modes lose their massless 
character, but only their "effective" number is reduced. Hence we believe that the spe- 
cific anisotropies we have treated here, besides their obvious importance for many real 
systems, supply a qualitative picture for what may happen rather generally. 

In order to show that by means of our renormalization procedure a coherent de- 
scription of the crossover from a critical theory to an asymptotically Gaussian model 
may be obtained, we discussed the time-dependent Ising model (or rather its field- 
theoretical representation) in the ordered phase. Although we have only pointed out 
the features of the single-component model for T < T c , it should be quite obvious how 
generalizations to other situations could be achieved. The effective couplings appearing 
in Wilson's flow functions vanish in the limit £ — > 0, and hence the power laws are 
eventually characterized by mean-field exponents. On the other hand, the expansion 
parameter of the perturbation series for the scaling functions themselves acquires a fi- 
nite, but non-universal value. Of course, none of the resulting integrals will contain any 
infrared divergence. This description remains valid until additional contributions come 
into play which are considered as irrelevant with respect to the critical properties and 
thus have been omitted in the Hamiltonian. 

After these preliminary discussions we have treated the complete crossover region 
for the relaxational models A and B with cubic anharmonicity. As for the isotropic sys- 
tem, the renormalization constants and flow functions were evaluated to one-loop order, 
and the ensuing flow equations were solved numerically. An important distinction has 
to be made between the cases n > n c = 4 (to one-loop order) , where a direct transition 
to the uncritical behavior takes place, and n < n c , where in an intermediate wavenum- 
ber and frequency range coexistence-type singularities appear. The latter crossover 
scenario crucially requires that in the vicinity of the critical point the 0(n)-symmetry 
is dynamically restored by fluctuation effects [25-27]. Lawrie's coexistence fixed point 
becomes unstable with respect to cubic distortions. This we have studied by using sev- 
eral values for the renormalized parameter y in Figs. 5, 7-9, although the latter is merely 
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a qualitative measure of the microscopic anisotropy. If y < 10~ 4 , all the characteristic 
features of Goldstone anomalies in three dimensions are still present in the effective 
static exponent 2 — rj^f (4.23), and they are even more prominent in the frequency 
dependence of the dynamic response and correlation functions. Therefore experiments 
investigating the low-frequency behavior of the dynamical susceptibility would be very 
enlightening: We emphasize that any hint of a minimum of the corresponding effective 
exponent (compare Figs. 7-9) would be a clear signature of the relevance of Goldstone 
modes. Table II provides an overview for the relevance of coexistence anomalies in the 
cubic case as the dimension d and the number of components n are varied. 

On the other hand, the dipol-dipol interaction induces an anisotropy in momentum 
space, which also explicitly breaks the 0(n)-symmetry. Yet the Goldstone modes are 
not entirely extinguished, but merely their number is in effect reduced by one. Hence 
while for n = 2 a crossover to an asymptotically uncritical theory takes place, for n > 3 
coexistence anomalies persist, governed by the dipolar fixed point (5.31a). Unfortu- 
nately the considerably more complex structure of the dipolar propagators (5.12) of 
model A renders a complete crossover theory a rather cumbersome task. However, the 
asymptotic model, characterized by diverging mass m(£) — > oo and dipolar coupling 
g(£) — > oo, may be considered on a very similar basis as in the isotropic case. Again, 
a one-loop theory for the two-point cumulants and the corresponding geometric series 
for the two-point vertex functions provides an exact representation in the ordered phase 
for q — > and u> — > 0. As a consequence of the anomalous dimension (5.31b) remaining 
unaltered, the asymptotic power laws for n > 3 are identical to those we have found for 
(7 = 0, see Table I. The anisotropy appears in the ^-dependence of the relevant scaling 
functions (5.37), in particular with respect to the angle between order parameter and 
external wavevector. If the vector Q of Eqs. (5.35b) or (5.38a,b), respectively, is not 
equal to zero, the amplitudes are generally smaller than in the isotropic case [compare 
(5.33)]. A qualitative summary of the different scenarios is given in Table III. 

We finally remark that all our results for the models discussed in this paper are 
in accord with exact theorems, especially concerning the possibility of long-range order 
in two dimensions [22,23]. Our formalism also allows extensions to more complicated 
models containing mode-coupling vertices, which usually constitute a relevant part of 
the dynamics of 0(n)-symmetric systems. In this paper we rather intended to study the 
influence of typical anisotropies than modifications by these reversible forces. Most of 
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the qualitative conclusions based on the relaxational models are expected to hold also for 
models with mode-coupling terms. For example, in ferromagnets the dipolar coupling 
g(£) diverges faster than the longitudinal mass parameter m(£) [compare Eq.(5.30)]. 
Hence near T c the influence of the dipolar interaction should dominate the effects due 
to the formation of a finite order parameter, as has been anticipated by Kotzler, Kaldis, 
Kamleiter and Weber in the interpretation of their experiments on EuS below T c [41]. 
A detailed analysis, however, requires a treatment of model J [1] appropriate for the 
dynamics of isotropic ferromagnets. Concerning the statics of such dipolar ferromagnets, 
we emphasize once again that the present theory is already complete in the coexistence 
limit. 
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APPENDIX: CUBIC TWO-POINT VERTEX FUNCTIONS 



In this Appendix we list the zero- and one-loop-diagrams (Fig. 10), and the corre- 
sponding analytical results for the two-point vertex functions of the relaxational models 
A (a = 0) and B (a = 2) including weak cubic anisotropy. The integration over the inter- 
nal frequency has already been performed via the residue theorem. (Compare Appendix 
B of Ref.[ll].) 

(-) = i«+W ^ W+ 2 -^ (mg-mg) J k ( ^ +fc2 / K+fc2) ) 

u 2 ^ (l + ^) a 1 

(6) = - A «• ^ ml J, m ,; (| _^ r ^ + (M - L7H g + ( | +£) 'j +( |_£ ) . Lmg + ( |_ £) -j 

(c) = -Ao ,« ^ m j fc _ 2+(| _ fc y ^ + (| + £)-[ m g + (| + £)-j + 

r aa(9,^) : 

(d) = ia>+Aog a (m 2 +q 2 ) 



2 f (f + fc)" i 

(/) - -Ao 3 6 ™ Jk rag + (| _g )a ^ + ( | +S) - [ m g + ($+£) a ] +(f-S)« K + d-A?) 2 ] 



roir7r(9,w) : 
(g) = -2A g a 

= 4 ™S Re j k H + {i+k y^ + [i _ m 



{i) = -2\ q a 

1 



(j) = -2 A q 2a ^ ^2. m 2 R e f, n . ^ ^, 2 - T i . ^ ^, 2 l X 



1 



+ +fc) a [rn 2 + (f+fc) 2 ] +(|-fc) a [^ + (|-fc) 2 ] 



^ + (|+fe) a [mg + (|+fc) 2 J +(|-fc)° [ml + {i-k) 2 \ 
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TABLE I. Coexistence anomalies of the isotropic relaxational models. 



Model A Model B 



,— e 



XL(q,0) ocq e oc q~ 

RexL(0,u/q a ) oc u~ e / 2 oc(u/q 2 )~ e / 2 
G L (0,u/q a ) ocu' 1 -"/ 2 (x(u/q 2 )~ e 



TABLE II. The influence of cubic anisotropy. 

d=2 d=3 d=4 

n = 1 crossover to a Gaussian theory, as for v = 
2 < n < 4 via — n —\ — *■ anomalies log. corr. 
n > 4 directly to 0, no anomalies log. corr. 



TABLE III. The influence of dipolar interactions. 

d=2 d=3 d=4 

n = 1 crossover to a Gaussian theory, as for g = 
n = 2 crossover to a Gaussian theory, no anomalies 
n > 3 = ^2 — *■ anomalies log. corr. 
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FIGURE CAPTIONS 



FIG.l. Flow of the effective coupling u e ff(£) for n = 1 and e = 1: (a) flow diagram, (b) 
universal crossover. 

FIG. 2. Effective exponent 2 — i] e ff of the static susceptibility in the case n = 1 and e = 1. 
FIG. 3. Basic elements of the dynamical perturbation theory for the cubic relaxational 
models below T c . 

FIG. 4. Schematic flow diagrams for m = (T = T c ) in the case n = 2 (a) and n = 8 (b); 
e = 1. 

FIG. 5. Flow of the effective coupling b e ff(£) for several values of the anisotropy parameter 
y (y = corresponds to the isotropic case); n = 2 and e = 1. 

FIG. 6. Schematic flow diagrams for m > (T < T c ) and v < in the case n = 2 (a) and 
n = 8 (b). 

FIG. 7. Effective exponents 2 — (a) and 2 — (b) of the static susceptibilities for 
several values of the anisotropy parameter y; n = 2 and e = 1. 

FIG. 8. Frequency dependence of RexL(0,o;) (a) and of Gl(0,i<;) (b) for the cubic model 
A with different values of the anisotropy parameter y; n = 2 and e = 1. 

FIG. 9. Frequency dependence of Re Xl(0, oj/q 2 ) (a) and of GL(0,Lv/q 2 ) (b) for the cubic 
model B with different values of the anisotropy parameter y; n = 2 and e = 1. 
FIG. 10. Zero- and one- loop diagrams for the two-point vertex functions of the cubic relax- 
ational models A (a = 0) and B (a = 2). The corresponding explicit analytical 
expressions are listed in the Appendix. 
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